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Second partial sums of m-th powers 

We consider the sums of powers of successive integers: 


n 

J2k m = l rn + 2 rn + ...+n m 

k = 1 

which, as we know, are calculated with the Faulhabcr polynomials, as fol- 
lows: 

n 

y^k — l/2(n 2 + n) 

k = 1 
n 

k 2 = l/6(2?i 3 + 3n 2 + n) 

k = 1 
n 

]T k 3 = l/4(n 4 + 2n 3 + n 2 ) 

k = 1 
n 

k 4 = 1/30 (6n 5 + 15n 4 + 10n 3 — n) ... etc. etc. 

k = 1 

Each of these formulas generates, as n varies, an increasing numerical 
sequence, of the type of that obtained for m — 2: 

1, 5, 14, 30, 55, 91, 140, 204, 285, ... 
that is the sequence of the square pyramidal numbers . 

We aim to find a way to calculate, given any of these sequences, the sum 
of its first n terms, that is, the second partial sums of the m-th powers. 

An opportunity to obtain this is offered by the following table: 
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By summing the content of each column (green+ycllow boxes), we obtain 
the sum of the first n m-th powers, that (in tribute to Faulhaber) we denote 
by F {m) : 


n 

F { m) = E k m 


k = 1 


2 





and the contents of the entire table will be then: 


( n + 1 )F( m ) — {n + 1) 

k= l 

The green section contains, in each row, the amount: K rn x K = K {m+V l 
By summing the contents of all rows of this section, one obtains the sum of 
the first n (m+l)-th powers: 

n 

F im+ 1) = E fc<”* +1) 

k = 1 

The yellow section contains, in the columns, the sequence of Fr m ) sums. By 
summing the contents of all columns of this section, one obtains the second 
partial sums of the m-th powers. 

The quantity that we seek is then obtained by subtracting to the content 
of the entire table, the content of the green boxes, that is: 


F{m) (jl “I” 1)-^(to) 


k = 1 


Ft 


(m+1) 


( 1 ) 


As we shall see later, this formula will enable us to conduct a research , 
limitless and in every direction, to find polynomial expressions generating 
integer sequences of sums of powers of every order and grade. 

By performing algebraic calculations for m = 3, you get: 

n 

J2 F( 3 ) = (n + 1 )F( 3 ) - ^( 4 ) = (n + l)[(n 2 + n)/2] 2 - (6n 5 + 15 n 4 + 

k = 1 

10n 3 - n)/30 = (3n 5 + 15n 4 + 25n 3 + 15n 2 + 2n)/60 
Polynomials obtained for m from 1 to 8 are listed in the following table: 


m 

Second partial sums of m -th powers 

i 

(n 3 + 3 n~ + 2 n )/6 

2 

(n 4 + 4n 3 + 5n 2 + 2n )/12 

3 

(3n 5 + 15 « 4 + 25 n 3 + 15« 2 + 2 n )/60 

4 

(2 n 6 + 12m 5 + 25 n 4 + 20« 3 + 3« 2 - 2 n )/60 

5 

(2 n 7 + 14 n 6 + 35 n 5 + 35n 4 + In 3 - In 2 - 2 n )/84 

6 

On 8 + 24n 1 + lOn 6 + 84« 5 + 21n 4 - 28 n 3 - 10« 2 + 4 n )/168 

7 

(5 n 9 + 45« 8 + 150n 7 + 210/7 6 + 63/7 5 - 105// 4 - 50/7 3 + 30/7 2 + 12/7 )/360 

8 

(2« 10 + 20n 9 + 75// 8 + 120n 7 + 42/7 6 - 84/! 5 - 50// 4 + 40/i 3 + 21/? 2 - 6 n )/180 
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Polynomial expressions generated by (1) are the natural extension of 
those listed at the beginning. The general formula for obtaining them in a 
direct way is written, in the compact notation of Faulhaber formula, in the 
following way: 


Y^F m =(n + l)F m -F (m+l) 

k = 1 

1 ™ . /m + 1\ 

7 X(-i)'( . K 

1 h=n & 


n + 

m + x k=0 


m+l—k 


1 , “+} /m + 2\ 

, >*(«+!) 

m + 2 


m+2—k 


k = 0 


where the quantities are the Bernoulli numbers . 

This formula calculates, for each natural number n, the second partial 
sums of the m-tli powers. 
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Partial sums of m-th powers with 
Faulhaber polynomials 

In previous section we saw how to use Faulhaber polynomials F( m ) to obtain 
formulas of second partial sums of m-th powers. 

We will seek now an iterative process for deriving, starting from the 
formula ( 1 ) of previous section, the polynomial expressions that calculate 
the partial sums later. 

Whereas each partial sum represents, by definition, the sum of the first 
q terms of the previous sum, we have, in general: 

q 

Pj(m) Pj-l(m) 

n= 1 

Then, the polynomial expression of the j-th partial sum of the m-th power 
is obtained by adding (from 1 to q) all terms of the polynomial of the previ- 
ous partial sum. In the practice is better to use non factored polynomials, 
ranked according to n powers: to each power of n will correspond the F( m ) 
polynomial. 

We will perform the derivation of the polynomial P 3 on the first power: 

P 2 ( 1 ) = (n 3 + 3 n 2 + 2n) / 6 

^ 3 ( 1 ) — 5Z -^ 2 ( 1 ) = ^(n 3 + 3 n 2 + 2n)/6 

n= 1 n= 1 

= (P(3) d" 3P(2) + 2F(!))/6 = n(l + n)(2 + n )( 3 + n )/ 24 

The knowledge^] of the general formula (1) enables then, using Faulhaber 
polynomials, the direct derivation of the m-th power formulas. 

With Mathematica , using mainly Factor and Table commands and a 
bit of copy and paste, you get faster: 


m 

Third partial sums of m-th powers: P 3 ( m j 

1 

n*(l + n)*(2 + n)*(3 + n)/24 

2 

n*(l + n)*(2 + n)*(3 + n)*(3 + 2*n)/120 

3 

n*(l + n)*(2 + n)*(3 + n)*(l + 3*n + n A 2)/120 

4 

n*(l + n)*(2 + n)*(3 + n)*(3 + 2*n)*(-l + 6*n +2* n A 2)/840 

5 

n*(l + n)*(2 + n)*(3 + n)*(-l + 2*n + n A 2)*(2 + 4*n + n A 2)/336 

6 

n*(l + n)*(2 + n)*(3 + n)*(3 + 2*n)*(2 - 30*n + 35*n A 2 + 30*n A 3 + 5*n A 4)/5040 

7 

n*(l + n)*(2 + n)*(3 + n)*(6 - 6*n - 20*n A 2 + 15*n A 3 + 25*n A 4 + 9*n A 5 + n A 6)/720 

8 

n*(l + n)*(2 + n)*(3 + n)*(3 + 2*n)*(l + 36*n - 69*n A 2 + 45*n A 4 + 18*n A 5 + 2*n A 6)/3960 




1 Ignoring the (1) we would not have P‘ 2 (\) that enabled us to derive P 3 (i). 
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m 

Fourth partial sums of m-th powers: P 4(m) 

i 

n*(l + n)*(2 + n)*(3 + n)*(4 + n)/120 

2 

n*(l + n)*(2 + n) A 2*(3 + n)*(4 + n)/360 

3 

n*(l + n)*(2 + n)*(3 + n)*(4 + n)*(2 + 4*n + n A 2)/840 

4 

n*(l + n)*(2 + n) A 2*(3 + n)*(4 + n)*(-l + 12*n + 3*n A 2)/5040 

5 

n*(l + n)*(2 + n)*(3 + n)*(4 + n)*(-24 + 20*n + 85*n A 2 + 40*n A 3 + 5*n A 4)/15120 

6 

n*(l + n)*(2 + n) A 2*(3 + n)*(4 + n)*(-l - 8*n + 14*n A 2 + 8*n A 3 + n A 4)/5040 

7 

n*(l + n)*(2 + n)*(3 + n)*(4 + n)*(48 - 100*n - 89*n A 2 + 160*n A 3 + 140*n A 4 + 36*n A 5 + 3*n A 6)/23760 

8 

n*(l + n)*(2 + n) A 2*(3 + n)*(4 + n)*(l + 4*n + n A 2)*(21 - 48*n + 20*n A 2 + 16*n A 3 + 2*n A 4)/23760 




m 

Fifth partial sums of m-th powers: P 5 ( m ) 

i 

n*(l + n)*(2 + n)*(3 + n)*(4 + n)*(5 + n)/720 

2 

n*(l + n)*(2 + n)*(3 + n)*(4 + n)*(5 + n)*(5 + 2*n)/5040 

3 

n*(l + n)*(2 + n)*(3 + n)*(4 + n)*(5 + n)*(10 + 15*n + 3*n A 2)/20160 

4 

n A 2*(l + n)*(2 + n)*(3 + n)*(4 + n)*(5 + n) A 2*(5 + 2*n)/30240 

5 

n*(l + n)*(2 + n)*(3 + n)*(4 + n)*(5 + n)*(-2 + 5*n + n A 2)*(9 + 10*n + 2*n A 2)/60480 

6 

n*(l + n)*(2 + n)*(3 + n)*(4 + n)*(5 + n)*(5 + 2*n)*(-3 + 5*n + n A 2)*(4 + 15*n + 3*n A 2)/332640 

7 

n*(l + n)*(2 + n)*(3 + n)*(4 + n)*(5 + n)*(-3 + 5*n + n A 2)*(-2 + 5*n + n A 2)*(5 + 5*n + n A 2)/95040 




... etc. etc. 
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Pascal’s triangle and recurrence relations for 
partial sums of m-th powers 

We build with Excel the table that calculates, for successive additions (left 
ccll+top cell) the partial sums of powers of positive integers: 
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2 th sums 
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i 
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We want to obtain the recurrence relation for second sums, that is a 
formula for calculating the n-th term in the column "2th sums" as a function 
of the previous terms. 

The formula that we seek is obtained by analyzing the data in the table 
as follows: 

c = a + b 
e = b + d 

f=c+e=a+b+e=a+e-d+e 
f = 2e - d + a 

Indicating with a(„ im ) the n-th term of the sequence, we therefore have: 

2th sums: a (n>m) = 2a (n _i, m ) - a (n _ 2 , m ) + n m 

Extending the previous scheme to the successive columns, one obtains: 

3th sums. ci(n,m) 3a (n _ 2im ) T ®(n— 3 ,m) T n 

4tll SUmS. U(n,m) 4a( n _i, m ) Q&(n—2,m) T U(n— 4 m ) T 71 

From this point forward we continue (successfully) using Pascal’s triangle , 
by alternating signs, with the following results: 

5tll SUmS. 5ct( n _i im ) 2 ,m) T 10ct(n— 3 ,m) 5(2( n _4 ^ m ) 

5,m) T Tl 

and so on ... 
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So, if we denote by p the order number of the partial sums, its recurrence 
relation is obtained by the following general formula: 


P_1 ^ 


a 


(n,m) 


K-i) 


k = 0 


ylc+lj 


a 


( n-l-k,m ) 


■ 1 
+ n 


This was a "non linear" recurrence relation between the sequence terms. 
There are also "linear" recurrence relations, obtainable from Pascal’s trian- 
gle, having the following general formula: 


p+m - 1 


a 


In 


.„> = I (-1) 


k = 0 


f p+m\ 


V k+l J 


a, 


( n—\—k,m ) 


Example: for the "fourth partial sums of sixth powers" you get the linear 
relationship: 

Oj[n) lOu^—p 45cq n _2) T 120(Z( ri _3) 210n^_4^ T 252(2^—5) 210tt( ri _(3) 

+120a( n _7) 45a( n _s) T 10ci( n _ 9 ) ffin— io) 

The coefficient lists of such expressions are shown in OEIS with the term 
"signature". 



Horizontal sequences 


We insert in Excel sequences of "Third partial sums of m-th powers", ar- 
ranging them in a table as follows: 


Second partial sums of m-th powers 

n 

o 

II 

§ 

m= 1 

m =2 

m =3 

m= 4 

m =5 

m— 6 

m —1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

2 

3 

4 

6 

10 

18 

34 

66 

130 

3 

6 

10 

20 

46 

116 

310 

860 

2446 

4 

10 

20 

50 

146 

470 

1610 

5750 

21146 

5 

15 

35 

105 

371 

1449 

6035 

26265 

117971 

6 

21 

56 

196 

812 

3724 

18236 

93436 

494732 

7 

28 

84 

336 

1596 

8400 

47244 

278256 

1695036 

8 

36 

120 

540 

2892 

17172 

109020 

725220 

4992492 

9 

45 

165 

825 

4917 

32505 

229845 

1703625 

13072917 

10 

55 

220 

1210 

7942 

57838 

450670 

3682030 

31153342 

11 

66 

286 

1716 

12298 

97812 

832546 

7431996 

68720938 

12 

78 

364 

2366 

18382 

158522 

1463254 

14167946 

142120342 


Consider the recurrence relation for 2th sums seen above: 


n,m ) l,m) ®(n— 2, in) 4“ ^ 


We will use this relationship to derive formulas of sequences that appear in 
each row of the table. 

The opportunity to do this lies in the fact that the above relationship 
is valid, for all n, in every column of the table , and then for each term of 
the considered horizontal sequence. You get formulas by replacing subse- 
quently, in the general recurrence relation, polynomials that gradually you 
get, starting from n = 2. This is an inductive process that applies to any 
order of partial sums of m-th powers. 

We perform here derivations of formulas until the 0(6, TO ) sequence. 


71 = 2 : 

71 = 3 : 

n=4: 

71=5: 

71=6: 


a (2 ,m) = 2x1-0 + 2 m = | 2 m + 2 


a (3 ,m) = 2 a (2 , m) - a (1 , m) + 3 m = 2(2 m + 2) - 1 + 3 m = 2 m+1 + 3 m + 3 


a ( 4,m) = 2a (3 , m) - a (2 , m) + 4 m = 2(2 m+1 + 3 m + 3) - (2 m + 2) + 4 m 



This process works indefinitely, generating polynomial expressions longer 
and longer, which in turn generate sequences with terms that magnify more 
and more rapidly. 1 personally found sequences up to the fifth partial sums, 
visible in the OEIS Encyclopedia to the above link. 
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Antidiagonal sums of j-th partial sums of m-th 
powers 

Let us examine the array formed by the j-th partial sums of third powers: 


j-th partial sums of third powers 

n 

1=0 

1=1 

1=2 

1=3 

1=4 

j=5 

j=6 

1=7 

i 

1 

1 

1 

1 

1 

1 

1 

1 

2 

8 

9 

10 

11 

12 

13 

14 

15 

3 

27 

36 

46 

57 

69 

82 

96 

111 

4 

64 

100 

146 

203 

272 

354 

450 

561 

5 

125 

225 

371 

574 

846 

1200 

1650 

2211 

6 

216 

441 

812 

1386 

2232 

3432 

5082 

7293 

7 

343 

784 

1596 

2982 

5214 

8646 

13728 

21021 

8 

512 

1296 

2892 

5874 

11088 

19734 

33462 

54483 

9 

729 

2025 

4917 

10791 

21879 

41613 

75075 

129558 

10 

1000 

3025 

7942 

18733 

40612 

82225 

157300 

286858 

11 

1331 

4356 

12298 

31031 

71643 

153868 

311168 

598026 

12 

1728 

6084 

18382 

49413 

121056 

274924 

586092 

1184118 


1 9 37 111 283 657 1441 3051 


An interesting sequence in this table is that obtained by summing data 
diagonally, as shown in figure, starting from the top left. You get: 1, 9, 37, 
111, 283, ... Inserting this sequence in "Mathematica”, using the command: 

F act or [Find Sequence Function[l, 9, 37, 111, 283, ..., n]] 

you get the formula: 

a(n) = 13(2” - 1) - 3n 2 - 9 n 

Inserting also sequences thus obtained from preceding and successive m-th 
powers tables, one obtains: 


m 


Formulas 

i 


1 

* (2 A n - 1 ) 


2 


3 

* (2 A n - 1) - 


3 


13 

* (2 A n - 1 ) - 

3 * n A 2 - 9 * n 

4 


75 

* (2 A n - 1 ) - 

4 * n A 3 - 18 * n A 2 - 52 * n 

5 


541 

* (2 A n - 1 ) - 

5 * n A 4 - 30 * n A 3 - 130 * n A 2 - 375 * n 

6 


4683 

* (2 A n - 1 ) - 

6 * n A 5 - 45 * n A 4 - 260 * n A 3 - 1125 * n A 2 - 3246 * n 

7 


47293 

* (2 A n - 1 ) - 

7 * n A 6 - 63 * n A 5 - 455 * n A 4 - 2625 * n A 3 - 11361 * n A 2 - 32781 * n 







Looking in OEIS for recurrences in this table, I discover that the blue 
sequence is given by A000670 . and A208744 gives the triangle of coefficients. 
This result speaks for itself, I plugged it at A257448 . 
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The differences of m-th powers 

The first differences of the m-th powers are computed, by definition, in the 
following way: 

Di(rn) = n m - (n- l) m (2) 

Using the program "Mathematica" you get from the (2) the following 
complete polynomials: 


in 

First differences of m-th powers: D 1(m) = n m - (n - 1)'” 


i 

1 


2 

-1 + 2 n 

OFF. 1 

3 

1 - 3 n + 3 n A 2 


4 

-1 + 4 n - 6 n A 2 + 4 n A 3 


5 

1 - 5n + 10n A 2 - 10n A 3 + 5n A 4 


6 

-1 + 6 n - 15 n A 2 + 20 n A 3 - 15 n A 4 + 6 n A 5 


7 

1 - 7 n + 21 n A 2 - 35 n A 3 + 35 n A 4 - 21 n A 5 + 7 n A 6 


8 

-1 + 8 n - 28 n A 2 + 56 n A 3 - 70 n A 4 + 56 n A 5 - 28 n A 6 + 8 n A 7 


9 

1 - 9 n + 36 n A 2 - 84 n A 3 + 126 n A 4 - 126 n A 5 + 84 n A 6 - 36 n A 7 + 9 n A 8 



You will get the subsequent differences by replacing n with (n — 1) in 
the corresponding polynomials having degree >1 of the preceding tables: 


m 

Second differences of m -th powers: D 2(m ) 


i 



2 

2 

OFF. 2 

3 

6 (-1 +n) 


4 

2 (7 - 12 n + 6 n A 2) 


5 

10 (-3 + 7 n - 6 n A 2 + 2 n A 3) 


6 

2 (3 1 - 90 n + 105 n A 2 - 60 n A 3 + 15 n A 4) 


7 

14 (-1 + n) (9 - 22 n + 23 n A 2 - 12 n A 3 + 3 n A 4) 


8 

2 (127 - 504 n + 868 n A 2 - 840 n A 3 + 490 n A 4 - 168 n A 5 + 28 n A 6) 


9 

6 (-1 + n) (85 - 296 n + 460 n A 2 - 408 n A 3 + 222 n A 4 - 72 n A 5 + 12 n A 6) 



m 

Third differences of m -th powers: D 3(m } 


i 



2 


OFF. 3 

3 

6 


4 

12 (-3 + 2 n) 


5 

30 (5 - 6 n + 2 n A 2) 


6 

60 (-3 + 2 n) (3 - 3 n + n A 2) 


7 

42 (43 - 90 n + 75 n A 2 - 30 n A 3 + 5 n A 4) 


8 

84 (-3 + 2 n) (23 - 42 n + 32 n A 2 - 12 n A 3 + 2 n A 4) 


9 

6 (3025 - 8694 n + 10836 n A 2 - 7560 n A 3 + 3150 n A 4 - 756 n A 5 + 84 n A 6) 
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m 

Fourth differences of m -th powers: D 4(m } 


i 



2 


OFF. 4 

3 



4 

24 


5 

120 (-2 + n) 


6 

120(13- 12 n + 3 n A 2) 


7 

840 (-2 + n) (5 - 4 n + n A 2) 


8 

168 (243 - 400 n + 260 n A 2 - 80 n A 3 + 10 n A 4) 


9 

504 (-2 + n) (185 - 272 n + 164 n A 2 - 48 n A 3 + 6 n A 4) 



and we could continue indefinitely . 

The so obtained polynomials, properly factored, provides formulas to 
calculate, at various orders, the differences of m-th powers. 

Note that, from each table, you get sequences whose offset is equal to 
the differences order number. 
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Conclusion 


In January 2015, I discovered that many of sequences obtained by the above 
procedures were not yet in the OEIS Encyclopedia database. This, in some 
way, gave me confirmation of newness of the introduced derivation proce- 
dures. I then immediately worked, contributing to the database OEIS with 
84 new integer sequences, all approved and published in about two months. 
You can consult them here . 
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